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(57) Abstract 

Methods and apparatus for delivery of ultrasound to the brain, without requiring removal of portions of the skull, call for transmission 
of ultrasound with a plurality of transducers (12) aimed to induce cavitation at least at a selected region of the brain. An excitation source 
(22) is arranged for driving at least selected transducers at differing phases with respect to one another, e.g., to compensate for phase shifts 
(or phase distortions) effected by the skull on the ultrasound output by each transducer. As a result, the ultrasound waves reaching the 
selected region from the transducers arrive substantially in phase with one another, e.g., within 90 degrees and preferably 45 degrees, and 
still more preferably 20 degrees of one another. 
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Methods and Apparatus for Delivery of Noninvasive Ultrasound 
Brain Therapy Through Intact Skull 



Sponsorship 



The research resulted, at least, in part, from work performed under NCI Research Grant 
No 46627. 

5 

Reference to Related Applications 

This application claims the benefit of the filing date of and is a continuation-in-part of 
copending, commonly assigned United States Patent Application Serial No. 08/71 1,289, filed 
10 August 21, 1996 (Attorney Docket No. 0092664-0008), the teachings of which are incorporated 
herein by reference. This application also claims the benefit of the filing date of copending 
United States Provisional Application Serial Nos. 60/034,084 (filed 12/23/96), and 60/045,453 
(filed 5/1/97). The teachings of those provisional applications are incorporated herein by 
reference. 
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Background of Invention 



The invention pertains to medical systems and, more particularly, to methods and 
apparatus for non-invasive application of focused ultrasound to the brain. The invention can be 
20 used, for example, in the diagnosis and treatment of neural ailments. 

According to the prior art, treatment of tissues lying at specific locations within the skull 
are limited to removal or ablation. While these treatments have proven effective for certain 
localized disorders, such as tumors, they involve delicate, time-consuming procedures that may 
25 result in destruction of otherwise healthy tissues. The treatments are generally not appropriate 
for disorders in which diseased tissue is integrated into healthy tissue, except in instances where 
destruction of the latter will not unduly effect neurologic function. 

The noninvasive nature of ultrasound surgery has special appeal in the brain where it is 
30 often desirable to destroy or treat deep tissue volumes without disturbing the healthy tissues. 
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Focussed ultrasound beams have been used for noninvasive surgery in many other parts of the 
body. Ultrasound penetrates well through soft tissues and, due to the short wavelengths (1.5 mm 
at 1 MHz), it can be focused to spots with dimensions of a few millimeters. By heating 
tumorous or cancerous tissue in the abdomen, for example, it is possible to ablate the diseased 
5 portions without significant damage to surrounding healthy tissue. 

Notwithstanding the potential benefits of ultrasound diagnostics and therapy of the brain, 
it has been commonly accepted that ultrasound cannot be applied through the intact skull. Early 
experiments, for example, showed that ultrasound is strongly attenuated by bone, and that brain 
10 tissue damage close to the skull results from the high temperatures caused by the energy loss. 
Accordingly, all of the ultrasound brain treatments performed so far have required the skull bone 
to be removed prior to the sonication. This makes the procedure invasive, and expensive with an 
added risk of complications. 

15 The requirement of surgical removal of the skull has been the main obstacle that has 

prevented ultrasound therapy to be widely tested in the brain despite the possibility of its clear 
benefits compared to other techniques. Accordingly, an object of this invention is to provide 
improved medical methods and apparatus, for diagnosis and therapy of the brain. A more 
particular object of the invention is to provide improved methods and apparatus for application 

20 of ultrasound to the brain. 

A more particular object of the invention is to provide such methods and apparatus as do 
not require removal of portions of the skull, via craniectomy or other such procedures. 

25 Still another object of the invention is to provide such methods and apparatus as can be 

used to precisely target regions within the brain. 

Still yet another object of the invention is to provide such methods and apparatus as can 
be used to effect heating or other physiologic change at such precisely targeted regions, without 
30 effecting substantial change in the surrounding, or other, regions of the brain or skull. 
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Another object of the invention is to provide such methods and apparatus as can be 
utilized over a wide range of ultrasonic frequencies. 

Still another object of the invention is to provide such methods and apparatus as can be 
5 implemented utilizing conventional materials. 

Yet still another object of the invention is to provide such methods as can be 
implemented without excessive expense. 
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Summary of the Invention 

The foregoing and other objects are met by the invention, which provides in one aspect 
methods and apparatus for delivery of cavitating ultrasound to the brain, without requiring 
5 removal of portions of the skull. 

Thus, in one aspect, the invention provides an apparatus for delivering ultrasound, 
through intact skull, to the brain comprising a plurality of transducers and an excitation source 
for driving each to induce cavitation at least at a selected region of the brain. The excitation 
10 source is particularly arranged for driving at least selected transducers at differing phases with 
respect to one another, e.g., to compensate for phase shifts (or phase distortions) effected by the 
skull on the ultrasound output by each transducer. As a result, the ultrasound waves reaching the 
selected region from the transducers arrive substantially in phase with one another, e.g., within 
90° and, preferably, within 45° and, still more preferably, within 20° of one another. 

15 

The excitation source drives the transducers to deliver ultrasound to the selected region at 
a frequency ranging from 0.01 MHz to 10 MHz and, preferably, from 0.1 MHz to 2 MHz. 
Sonication duration for the ultrasound ranges, according to further aspects of the invention, from 
100 nanoseconds to 30 minutes. According to still further aspects, the invention provides for 
20 delivery of ultrasound to the selected region with continuous wave operation or burst mode 
operation, where burst mode repetition varies from 0.01 Hz to 1 MHz. 

Further aspects of the invention provide an apparatus as described above, in which only a 
single transducer is used. 

25 

Still further aspects of the invention provide methods for operating transducer arrays as 
described above. 

These and other aspects of the invention are evident in the drawings and in the text that 
30 follows. 
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Brief Description of the Drawings 

A further understanding of the invention may be attained by reference to the drawings, in 

which: 

5 

Figure 1 depicts an embodiment of the invention and an experimental setup for testing it. 

Figure 2 depicts an embodiment of the invention for application of ultrasound to the brain 
of an animal. 

10 

Figure 3 depicts a phased array for application of ultrasound to the brain in accord with one 
practice of the invention. 

Figures 4A-4H illustrate the ultrasound pressure amplitude distribution in water across the 
15 focus of a transducer according to the invention at various frequencies, with and without skull 
sections in front of the transducer. 

Figures 5A and 5B illustrate the effect of applying ultrasound in accordance with the 
invention to brain tissue. 

20 

Figures 6A and 6B illustrate phase errors measured at the focus of ultrasound transducer 
arrays with a piece of skull in front of the transducers. 

Figures 7A-7C illustrate the pressure amplitude profiles across the focus of an ultrasound 
25 transducer phased array in water, through the bone, and through the bone with phase correction. 

Figure 8 illustrates the pressure amplitude distribution along the central axis of an ultrasound 
transducer array without and with the phase correction. 
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Figures 9A-9C illustrate the ultrasound pressure amplitude distribution measured across the 
focus of an ultrasound phased array in water, through skull without phase correction, and through 
skull with phase correction. 

Figure 10 depicts an embodiment of the invention for delivery of cavitating ultrasound to a 
patient's brain through the skull using a multi-element transducer array. 

Figure 1 1 depicts a method for delivery of cavitating ultrasound to a patient's brain through 
the skull using a transducer array. 
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Detailed Description of the Illustrated Embodiment 

Discussed below are methods and apparatus according to the invention for noninvasive 
delivery of ultrasound through intact skull to the brain. These permit ultrasound propagation 
5 through skull to effect cavitation, without causing undesired heating of the brain or, more generally, 
central nervous system (CNS) tissues. These also deliver adequate ultrasound power to ablate 
tissues, or to otherwise induce changes, at focal points (or regions) within the brain. As used herein, 
"tissue" refers to fluids, tissues or other components on or within a patient's body. 

10 Figure 10 depicts an apparatus according to the invention for delivery of ultrasound to the 

brain. The apparatus 10 includes an array of transducers 12 disposed on or near the external surface 
of the head of a human patient. The array 12 can constitute a single transducer, e.g., a spherically 
curved piezoelectric bowl of the type described below, though preferably, array 12 comprises a 
plurality of transducers arranged in a one-, two- or three-dimensional configuration. 

15 

Referring to Figure 3, for example, in one embodiment of the invention, array 12 comprises 
60 individual piezoelectric ceramic transducers mounted in a bowl of circular cross-section. The 
transducer elements, which can be, for example, 1 cm 2 piezoelectric ceramic pieces, are mounted 
in silicone rubber or any other material suitable damping agent for minimizing the mechanical 

20 coupling therebetween. Transducer arrays of this type are known in the art, as described, for 
example in Fan et al, "Control of the Necrosed Tissue Volume During Noninvasive Ultrasound 
Surgery Using a 16-EIement Phased Array," Medical Phvsics . v. 22, pp. 297 et seq (1995), a copy 
of which is filed as an appendix hereto and the teachings of which (e.g., at Figure 1 and the 
accompanying text) are incorporated herein by reference. The construction of a spherically curved 

25 phased array comprising multiple square-element transducers is shown in Figure 1 and the 
accompanying text of that publication. 

In the illustrated embodiment, each transducer of array 12 is independently driven by power 
and the control elements 18-22 to generate ultrasound for transmission through the patient's skull 
30 into the CNS tissues. More particularly, the transducers in array 12 are individually coupled, via 
coaxial cables 16, to separate channels of a driving system 18. Each channel of that system 18 
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includes an amplifier and a phase shifter, as shown. A common radio frequency (RF) signal is 
driven to each channel by radio frequency generator 22. Together, the radio frequency generator 
22 and driving system 18 drive the individual transducers of array 12 at the same frequency, but at 
different phases, so as to transmit a focused ultrasound beam through the patient's skull to a selected 
5 region within the brain. Unlike prior art systems, there is no need to remove portions of the skull 
beneath the array 12, e.g., via craniectomy or other such surgical procedure. 

The radio frequency generator 22 can be of any commercially available type. A preferred 
such generator is available from Stanford Research Systems, Model DS345. The generator is 
10 operated in a conventional way so as to generate an excitation signal, which is amplified and phase- 
shifted by the individual channels of driving system 18, in order to induce the corresponding 
transducers of array 12 to radiate ultrasound (e.g., in the range 0.01 MHz to 10 MHz). 

As illustrated, each channel in the driving system 18 includes a radio frequency amplifier. 
15 These can be any RF amplifiers of the type commercially available in the art. 

The phase shifting component of each channel of driving system 18 is constructed and 
operated in the conventional manner known in the art. Particularly, each phase shifter shifts the 
phase of an incoming RF excitation signal, received from RF generator 22, by an amount a u a 2 , a 3 , 
20 etc., as shown in the drawing. These phase shift factors a b a 3 , etc., can be pre-stored in the 
channels of driving system 18 or, preferably, generated by a controller 20. That controller 20 can 
be a general purpose, or special purpose, digital data processor programmed in a conventional 
manner in order to generate and apply phase shift factors in accord with the teachings hereof. 

25 The phase shift factors, a„ a* a* etc. serve two purposes. The first is to steer the composite 

ultrasound beam generated by transducer array 12 so that it is focused on a desired region within the 
patient's brain. The component of each phase shift factor associated with steering is computed in 
the manner known in the art for steering phased arrays. See, for example, Buchanan et al y 
"Intracavitary Ultrasound Phased Array System," TEEE Transa ctions Biomedical Engineering, v. 

30 41, pp. 1 178-1 187, a copy of which is filed as an appendix hereto and the teachings of which are 
incorporated herein by reference. Array steering, or focusing, is particularly discussed in that article, 
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for example, at pages 1179-1181 and, more particularly, in the section entitled "Focusing 
Techniques," the teachings of which are incorporated herein by reference. 

The second component of each phase shift factor a u a* tx 3 , etc., compensates for phase 
5 distortion effected by the skull in the ultrasound ouput by each transducer. In other words, the 
second component of the phase shift factors compensates for perturbations and distortions 
introduced by the skull, the skin/skull interface, the dura matter/skull interface, and by variations 
in the skull thickness. As those skilled in the art will appreciate, the two components that make up 
the phase shift factor for each channel of the driving system 18 are summed in order to determine 
10 the composite phase shift factor for the respective channel. 

The phase corrections that constitute the aforementioned second component of each 
phase shift factor can be determined a number of ways. In one embodiment of the invention, that 
component is determined from measurements of the thickness of the patient's skull under each 
1 5 transducer in array 12. Such skull thickness measurements can be made using conventional 
imaging techniques, such as computed tomography (CT) or magnetic resonance imaging (MRI). 

In an alternative embodiment, the aforementioned second component of each phase shift 
factor is determined by placing the array 12 on the patient's head and exciting individual 
20 transducers with a short ultrasound pulse. The echo back from the inner surfaces of the skull are 
monitored by the transducer array 12. The effect of the skull on ultrasound generated by each 
transducer is determined from those echos in accord with conventionally known relations. 

In still further alternative embodiments, the aforementioned second component of each 
25 phase shift factor is determined by implanting small hydrophones in the patient's brain. These 
are used to monitor the phase of the ultrasound generated by each transducer, e.g., in a manner 
similar to that described below in connection with Figure 1 . 

30 In lieu of illustrated components 18-22, the transducer array 12 can be driven by a 

driving system of the type disclosed in Buchanan et al, supra e.g. at Figure 2 thereof, the 
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teachings of which are incorporated herein by reference. Such a driving system would, of 
course, require modification in accord with the teachings hereof in order to incorporate phase 
shift factors a u a 2 , a 3 , etc., having first and second components as described herein and above. 

5 Referring to Figure 1 1, the system 10 is operated as described below in order to deliver 

ultrasound through the patient's skull to induce cavitation at a desired region of the brain. 

In step 24, the transducer array 24 is positioned on the patient's head. This is preferably 
accomplished in the conventional manner known in the art for insuring ultrasound transmission 
10 to the brain. The array is typically positioned over, and as close to, the region in which 
cavitation is to be induced. However, where intervening or adjacent cranial or CNS tissues 
might be adversely affected, the array can be positioned elsewhere and focused accordingly. 

In step 26, the aforementioned second component of the phase shift factor for each 
15 transducer is determined. This is accomplished in the manner described above, e.g., by 
individual exciting each element of the array and measuring the echo back. The alternative 
mechanisms described above can also be used to determine those components. Those skilled in 
the art will appreciate that in instances where the alternative mechanisms are used, they need not 
be performed after the array is positioned but, can be performed at some other prior time. 

20 

In step 28, the remaining components of each transducers' phase shift factor are 
determined. Particularly, those components associated with steering the array for delivery of 
ultrasound to the desired region are determined. Such determination is made, as indicated above, 
in the conventional manner known in the art for steering phased arrays. 

25 

In step 30, the array is excited, e.g., by control and driving elements 18-22, to focus 
ultrasound in the patient's head. As noted throughout, because the invention provides correction 
for phased distortion induced by the skull, that ultrasound can be supplied directly through the 
skull without the need for removal of a piece thereof. The ultrasound is applied in doses and 
30 timing sufficient to induce cavitation in the desired region, which may be, e.g., from 1 mm 3 - 1 
cm 3 , or larger. Those skilled in the art will appreciate that ultrasound waves in the frequency 
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range of 0.01 MHz to 10 MHz and, preferably, from 0. 10 MHz to 2 MHz can be applied with 
sonication duration ranging from 100 nanoseconds to 30 minutes, with continuous wave or burst 
mode operation. The burst mode repetition varies from 0.01 Hz to 1 MHz. 

. 5 In embodiments where the transducer array 12 includes only a single transducer, e.g., a 

10 cm diameter piezoelectric ceramic element as described elsewhere herein, step 26 is not 
utilized. In such an embodiment, the "array" is aimed based on its focal point. This is 
determined as a function of the size, radius of curvature and frequency output of the transducer 
in the manner known in the art. In a preferred embodiment, these factors are adjusted so that the 

10 transducer can be placed directly on the patient's skull, as above. However, where minor 
corrections are necessary, the transducer can be spaced apart from the skull, as necessary, in 
order to insure proper positioning of the focal point. 

Theory 

15 When an ultrasound beam propagates to a deep target location in the brain part of the 

energy is reflected back at the skin-skull interface due to the high acoustic mismatch between 
these two tissues. The propagating wave in the skull suffers attenuation losses due to scattering 
and absorption. The acoustic mismatch at the bone-dura interface causes part of the remaining 
wave to reflect back to the skull. The total insertion loss through skull depends on the 

20 frequency, and can be, on average, about 10 dB, and 20 dB at 0.5, and 1.5 MHz, respectively. 
The wave is further attenuated by absorption (about 5 Np/m/MHz) while it travels through the 
brain to the target volume. 

An ultrasound beam delivered to the brain can effect change in CNS tissues and fluids 
25 (herein, simply "tissue" or "brain tissue," etc.) by two mechanisms: heating and cavitation. The 
ultrasound beam can heat the tissue temperature due to energy absorption from the wave 
resulting in different degrees of thermal damage to the tissue depending on the temperature 
reached. For exposures of a few seconds, temperatures of above about 60° C are adequate to 
coagulate proteins and thus, necrose the tissue. The induced temperature elevation during short 
30 ultrasound exposures depends mainly on the absorbed power (<q>) although the shape and size 
of the focal spot can have a significant impact due to thermal conduction. The rate of 
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temperature rise (dT/dt) at the very beginning of an ultrasound pulse can be calculated from the 
pressure amplitude of the field (P), as follows: 

dT/dt = <q>/pc, [1] 

5 

where 



<q>= ccPV pv [2] 

p is the density (p b = 1030 kg/m 5 , p, = 1380-1810 kg/m 3 ) (s= skull, b = brain), 
10 c is the specific heat of the medium (c, =3.9 kJ/kg/°C ; c, =2.1-2.7 kJ/kg/°C), 

v is the speed of sound, and 

a is the amplitude absorption coefficient of the tissue (a b = 5Np/m/MHz, a, = 50 
Np/m at 0.5 MHz, a 8 = 300 Np/m at 1.5 MHz (if all attenuated energy is assumed 
to be absorbed)). 

15 

In order to achieve the same temperature within the target volume as in the skull, the 
ultrasound beam has to be focused to overcome the difference in the acoustic properties. The 
square of the pressure amplitude (P 2 ) is directly proportional to the ultrasound beam area 
allowing the required area gain (AG) to be calculated from Equations [1] and [2] by making the 
20 rates of temperature rise equal in the skull and the brain: 

AG = (P b VP, 2 ) = (a/a b )(p b /p.y(v b /v.)(cyc l ) = 30 at 0.5 MHz [3] 

60 at 1.5 MHz 

25 In addition the area gain has to compensate for the energy loss due to the skull and 

attenuation in the brain between the skull and the focal point: 

AG- fa e- 2P& )- 1 = 13 at 0.5 MHz [4] 

250 at 1.5 MHz 

30 
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where, 

cti is the insertion loss of skull, 

P is the amplitude attenuation coefficient, 

f is the frequency and x is the depth in the brain. 

5 

The total area gain is the product of these two area gains and is approximately 400 and 
15000 at 0.5 MHz and 1.5 MHz, respectively when the focus is located at the depth of 6 cm in 
the brain. These calculations are first order estimates and do not take into account phase shifts 
introduced by the variable thickness of the bone or thermal conduction and perfusion effects. 

10 

The second mechanism, cavitation requires negative pressure amplitudes that are large 
enough to form gas bubbles in the tissue. The pressure wave causes the bubbles to expand and 
then collapse. The collapse of the bubbles causes high temperatures and pressures that can cause 
direct mechanical damage to the tissue. Cavitation can offer more therapeutic options than 
15 thermal exposures of brain. The cavitation threshold in the soft tissues and in bone appears to be 
similar. 

Thus, focusing a cavitation-inducing ultrasound beam overcomes the attenuation losses 
in the bone and brain, but need not overcome differences in absorption coefficients, as is the case 

20 with the heat-inducing exposures. The beam area of cavitation-inducing ultrasound propagating 
through the skull has to be about 13 and 250 times larger than the focal area at frequencies of 0.5 
and 1.5 MHz, respectively. These area gains are 30 and 60 times smaller than the gains required 
for induction of thermal effects. The cavitation is not influenced by thermal conduction or 
perfusion effects. Therefore, it is clear that cavitation has significant advantages over the 

25 thermal effects. This is particularly true in instances where the ultrasound energy must be 
delivered to small focal regions that require high frequencies. 

Cavitation requires high pressure amplitudes but only short exposure durations, therefore 
cavitational effects can be induced without significant temperature elevation. For CNS tissues, 
30 sonications with durations of only 1 ms are adequate for bubble formation. The required peak 
intensities at 0.936 MHz during these sonications are measured to be around 4000 Wcnr 2 and 
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2000 Wcrrr 2 at 1 ms and 1 s exposures, respectively. Using these intensity values and the 
average ultrasound attenuation in brain of 5 Np/m and in bone 120 Np/m at 1 MHz, the 
maximum peak temperature elevation in the brain can be estimated (from equations 1 and 2) to 
be about 60° C and 0.1° C during the 1 s and 1 ms exposures, respectively. The corresponding 

5 temperature elevations in bone, if the focus was in the bone, are 1800° C and 3.6° C. During the 
sonications the temperature elevation in the bone would be reduced proportionally with the area 
gain. These values are frequency dependent. For example, bone heating would be about 13° C 
for 1 ms pulse at 1.5 MHz. This short thermal exposure is below the threshold for tissue 
damage. Thermal exposures can be further reduced using multiple pulses that can be repeated at 

10 a low frequency (for example 0. 1 Hz) thus, eliminating a temperature build up. 

Examples 

In one exemplary embodiment of the invention, which is illustrated in Figure 1, 
15 ultrasound was generated using a one-element transducer array having a spherically curved 10 
cm diameter piezoelectric ceramic (PZT4) bowl mounted in a plastic holder using silicon rubber. 
The ceramic had silver or gold electrodes both on the front and back surface. The electrodes 
were attached to a coaxial cable that was connected to a LC matching network that matched the 
electrical impedance of the transducer and the cable to the RF amplifier output impedance of 50 
20 ohm and zero phase. The matching circuit was connected to an RF- amplifier (both ENI A240L 
and A500 were used in the tests). The RF signal was generated by a signal generator (Stanford 
Research Systems, Model DS345). 

The ultrasound pressure wave distributions were measured using needle hydrophones 
25 (spot diameter 0.5 and 1 mm) and an amplifier (Precision Acoustics Ltd). The amplified signal 
was measured and stored by a oscilloscope (Tektronix, model 243 1L ). The hydrophone was 
moved by stepper motors in three dimensions under computer control. The pressure amplitudes 
measured by the oscilloscope were stored by the computer for each location. 

30 A piece of human skull (top part of the head: front to back 18 cm and maximum width 12 

cm) was obtained and fixed in formaldehyde. The acoustic properties of formaldehyde fixed 
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skull and a fresh skull are almost identical. The ultrasound transducer under test was positioned 
in a water tank the walls and bottom of which were covered by rubber mats to reduce ultrasound 
reflections. The tank was filled with degassed deionized water. The hydrophone was connected 
to the scanning frame, and positioned at the focus of the ultrasound field. 

5 

Utilizing the setup shown in Figure 1, the embodiment was tested at four different 
ultrasound frequencies: 0.246 MHz, 0.559 MHz, 1 MHz, and 1.68 MHz. The maximum peak 
pressure amplitudes achievable through the skull at the focus of the transducer was measured at 
each frequency. A shock wave hydrophone (Sonic Technologies Inc, ) was positioned at the 
- 10 acoustic focus. Bursts of 10-20 cycles were used to separate the acoustic signal from the 

electrical interference that was picked up by the hydrophone during sonication. Results of the 
testing are shown in Figures 4 A-4H. 

0 Particularly, Figure 4A illustrates the ultrasound pressure amplitude distribution in water 
15 at the focal point of the single transducer driven at 0.246 MHz, without the skull section in place. 
Figure 4B illustrates this same distribution when the skull was positioned in front of the 
transducer as illustrated in Figure 1. Figures 4C and 4D illustrate the same distributions (i.e., 
with and without the skull section in place) for a frequency of 0.559 MHz. Figures 4E and 4F 
illustrate the same distributions for a frequency of 1 MHz. Figures 4G and 4H illustrate the same 
20 distributions for a frequency of 1 .68 MHz. 

To demonstrate that the single-element array of this embodiment delivers sufficient 
energy to induce tissue damage through the skull, in vivo rabbit experiments were performed. In 
these experiments a window of about 15x15 mm was created in the top of the skull. The skin 

25 was placed over the skull opening and the animal was allowed to recover. A minimum of two 
weeks after the surgery the animal was anesthetized again and placed on top of a sonication tank 
as illustrated in figure 2 and the 0.556 MHz transducer was aimed such than the focus was 
located at 10-15 mm in the brain. The skull piece was positioned in between the traducer and the 
animal. A thermocouple probe (0,05 mm constantan and copper wires were soldered together at 

30 the tip) was placed on the skull bone (on the side of the transducer that is expected to be the 
hottest location) under a thin layer of connective tissue that was still attached on the skull. Then 

15 
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10 sonications at the maximum power level for the duration of 0.2 s were repeated with the rate 
of 1 Hz. The animal position was moved and the sonication repeated four times in the same 
location with a delay of about 5 min between the sonications to allow the bone temperature to 
return to the baseline. During the 10 s of pulsed sonication the bone temperature increased from 
5 the baseline of about 30° C to maximum of 43° C with rapid decay. After the sonications the 
rabbit was taken to a MRI scanner and Tl, T2 and contrast enhanced scan were performed. 
After the imaging the animal sacrificed. 

Figure 5 A is a scan of the rabbit brain illustrating the effect of 10 sonications for the 
10 duration of 0.2 seconds, with a pressure amplitude of 8 MPa, repeated at a rate of 1. Hz. During 
the 10 seconds of pulse sonication, the bone temperature in the rabbit skull peaked from a 
baseline of about 30° C to a maximum of 43° C, with rapid decay. The Figure is a T2-weighted 
fast spin echo image across the brain. The arrow in the Figure shows tissue damage at the focal 
point of the transducer. The skull window on the top of the head is facing down and, thus, the 
15 ultrasound beam propagated from bottom up. Figure 4B is identical to Figure 4 A, except insofar 
as it shows the results where the above sonication was repeated four times. 

Further embodiments of the invention utilize multi-element phase arrays of the types 
illustrated in Figures 3 and 10, in lieu of a single transducer. By controlling the phase of the 
20 ultrasound wave as a function of transducer location, these embodiments eliminate the phase 
distortion caused by the skull and thus, allow accurate aiming and use of higher frequencies, 
thus, permitting application of ultrasound to induce cavitation through the intact skull in regions 
of 1 mm to 1 cm 3 . 

25 Two phased arrays comprising these further embodiments hkd similar structure and the 

same driving hardware; the resonant frequency being their only significant difference. The two 
arrays operated at 0.6 MHz and 1.58 MHz. The radius of curvature of both of the transducers 
was 10 cm and both of them were cut into approximately 1 cm 2 square elements, as shown in 
Figure 3. The total number of elements in both arrays was 64 although only 60 were driven in 

30 the experiments due to hardware limitations. The ceramic bowl was cut using a diamond wire 
saw so that the elements were completely separated by a 0.3-0.5 mm kerf. The kerf was filled 
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with silicone rubber that kept the array elements together and isolated them acoustically. The 
silicone rubber allowed the transducer elements to vibrate with minimum amount of clamping. 
Each transducer element was connected to a coaxial cable and a matching circuit that was 
individually tuned. The arrays were similar to the one described in Fan et al, supra, at Figure 1 
5 and the accompanying text, the teachings of which are incorporated herein by reference. The 
array was driven by an in-house manufactured 64 channel driving system that included an RF 
amplifier and phase shifter for each channel. The phase and amplitude of the driving signal of 
each channel was under computer control, as described in Buchanan et al, supra, e.g., at Figure 2 
and the accompanying text, the teachings of which are incorporated herein by reference. 

10 

In addition to phased arrays configured as described herein and shown in the 
accompanying drawings, phased arrays can also be constructed in accord with the arrangements 
described and shown in co-pending, commonly assigned patent application 08/747,033, filed 
November 8, 1996, the teachings of which are incorporated herein by reference. 

15 

Using the methodologies and apparatus described above, it is possible to produce a sharp 
focus through the skull with the single element transducer when the operation frequency was 1 
MHz or less. The beams had secondary peaks introduced by the skiill but the main peak was the 
largest. The location of the peak was shifted by the skull by 1-3 mm from the geometric focus, as 
20 shown in Figures 4A - 4H. However, the focus was completely obliterated with an operating 
frequency of 1 . 67 MHz, as shown in Figure 4H. 

The maximum pressure amplitudes achieved at the maximum power output of the 
amplifier were frequency dependent and are given in the table below. The maximum average 
25 pressure amplitude at the frequency of 0.554 MHz was 8.0+/- 0.6 MPa. 
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10 



Average Pressure P+ P- 
Frequency (MHz) Amplitude (MPa^ <MPa) (MPa) 

0.248 3.8 4.2 3.3 

0.559 8.0 8.9 7.1 

1.0 5.9 6.9 4.5 

Figure 5 shows the image across the brain for the first of the sonications and demonstrate 
tissue damage indicated T2 changes. The tissue damage was also visible in Tl images with and 
without contrast enhancement. 



15 To measure the phase distortion caused by the skull, a hydrophone was placed in the 

geometric focus of the array under test. The skull was placed between the array and hydrophone 
and each transducer element was powered separately in sequence while recording the time 
difference between the reference signal and the acoustic wave at the focus. This was done with 
both of the arrays. The phase changes required to correct all of the waves to arrive at the same 

20 phase at the focus are plotted in Figure 6. 

To investigate the effect of the phase correction the pressure amplitude distributions were 
measured in water by scanning the needle hydrophone. The main impact of the phase at 0,6 MHz 
was in the location of the focus which could be corrected back to the geometric focus. This is 

25 shown in Figure 7A, which illustrates the pressure amplitude profile across the focus of the 0.6 
MHz phased array in water. Figure 7B shows the pressure amplitude profile across the focus 
through bone. And, Figure 7C shows the pressure amplitude profile through bone when a phase 
correction according to the invention is used. Figure 8 likewise illustrates the pressure amplitude 
distribution along the central axis of the array with and without phased correction. The 

30 magnitude was reduced to 33 % and 40 % of its water value without and with the phase 

correction, respectively. The importance of the phase correction was demonstrated more clearly 
with the higher frequency array. With this array the focus was completely destroyed by the skull 
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(Figure 9b). However, when phase correction was introduced, the focal spot was returned into 
its original shape (Figure 9c) with the half-width of the focus of about 1 mm. The insertion of the 
skull reduced the peak pressure at the focus of the 1.58 MHz array to about 5 % of its water 
value when phase correction was applied. 

5 . 

The embodiments of the invention discussed above and shown in the drawings provide 
improved methods and apparatus for neural diagnosis and therapy through application of short, 
high intensity ultrasound beams that induce cavitation at selected locations within the brain. 
These and other embodiments can be beneficially used to deliver focused ultrasound beams to 

10 the CNS tissues and fluids, thereby, permitting their ablation or other physiological modification. 
Thus, for example, the embodiments can be used to ablate tumors, cancers and other undesirable 
tissues in the brain. They can also be used, for example, in connection with the technologies 
disclosed in copending, commonly-assigned U.S. Patent Application No. 08/71 1,289 (the 
teachings of which are incorporated herein by reference) for modification of the blood-brain 

15 barrier, e.g., to introduce therapeutic compounds into the brain. Because they do not require that 
portions of the skull be removed, the embodiments permit the foregoing to be performed 
noninvasively. 

The results also show that adequate ultrasound transmission can be induced through 
20 human skull to induce cavitation in vivo. This can be done with single element applicators, e.g., 
preferably at frequencies less than 1 MHz and at higher frequencies with phased arrays that 
correct the phase distortion caused by the variable thickness of the skull. 

The maximum pressure amplitude of 8 MPa induced through the skull at 0.559 MHz was 
25 able to induce cavitation damage in vivo rabbit brain. This value was reached through an area of 
10 cm in diameter to a focal spot diameter of about 5 mm (50 % beam diameter). If the whole 
available skull surface around the brain is utilized, then a window of at least three times larger 
could be used. In addition, the geometric gain would allow the peak power through the skull to 
be increased. Acoustic power up to 30- 80 W/cm 2 of the transducer surface area for continuous 
30 wave sonication can be generated by ceramic transducers. Higher peak powers could be achieved 
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with the pulsed sonication used for induction of cavitation. Thus, it is estimated that much higher 
pressure amplitudes than measured here can be induced in the brain through the skull. 

The values measured in connection with the foregoing compare favorably with the 4 MPa 
5 that was reported to be the threshold value in vivo muscle at 0.6 MHz and a value of 8.5 MPa at 
0.936 MHz in vivo rabbit brain (the threshold at 0.6 MHz would be lower since it has been 
shown to decrease with frequency). Thus, the results demonstrate that adequate ultrasound 
transmission through skull can be generated to induce cavitation in the brain. 

10 Our results demonstrate that low frequency beams can be focused through the skull, 

though, the focus may be shifted from its geometric location. Therefore, it can be helpful to 
detect the focal spot location in the brain prior to the therapy exposure, e.g., using magnetic 
resonance imaging to detect the local temperature elevation or cavitation in the brain at exposure 
conditions that are below the tissue damage threshold. For example, low power test exposures 

15 can be delivered through the skull while using MRI to detect the location of the focal spot. 
Based on the imaging information the location can be corrected to overlap the target volume 
prior to the therapeutic exposure. 

The phase measurements with the arrays support the observation made with the single 
20 element transducers showing that at 0.6 MHz 80 % of the phase errors caused by skull are less 
than 90° and thus, each wave is adding to the pressure wave at the focus. However, at 1.58 MHz 
over half of the waves had phase shifts that caused the waves to arrive out of phase at the focus. 
This observation can be explained by the difference in wavelength that is 2.50 mm at 0.6 MHz 
and 0.95 mm at 1.58 MHz. 

25 

The possibility of inducing selective thermal damage at the focus, without damaging the 
skin or brain surface, may be possible due to the small focal spots achieved with the phase 
correction. However, the thermal exposures have to be short to reduce blood flow and perfusion 
effects that are strong in brain tissue. The sharp temperature gradients at the focus transport 



20 



WO 98/07373 



PCT/US97/14760 



more energy away from the focus than in the bone where the beam is wide and the gradients 
shallow. At 1 .58 MHz, full utilization of the skull surface may provide marginally adequate 
geometric gains to overcome the skull heating problem. However, at lower frequencies 
especially around 0.5 MHz the focal brain tissue thermal therapy seems feasible although not as 
5 likely as utilization of cavitation effects. 

The results demonstrate that the effects of the skull to the beam shape can be eliminated 
using a phased array with proper phase corrections. In the example above, the phase correction 
was calculated from hydrophone measurements. The same corrections can be made by 

10 measuring the skull thickness from CT or MRI scans and then calculating the phase correction 
required for each array element. The same may be accomplished by sending a short ultrasound 
pulse from each or selected elements of the of the phased array and then listening for the echo 
back from the inner surfaces of the skull or other structures in the brain. The effect of the skull 
on the wave propagation at each location could then be calculated. This can also be done before 

15 therapy by mapping the skull effect using ultrasound. 

Although good results were achieved with only 60 transducer elements in the phased 
array, still more and smaller elements may facilitate moving the focal spot inside of the brain. 
Similarly, at higher frequencies, smaller elements may allow better phase correction further 
20 reducing the losses induced by the skull. 

Thus, the geometric gain of about 20 that is required to compensate for the losses caused 
by the skull can be easily achieved by focusing. This is larger than the gain of 10 required to 
compensate the average losses. This indicates that adequate power for induction of cavitation 
25 can be delivered using phased arrays through the skull even at frequencies that are too high with 
a single element applicator. 

In summary, the invention provides methods and apparatus for noninvasive diagnosis and 
treatment of the brain using cavitational mechanism and pulsed ultrasound. It permits adequate 
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power transmission through the human skull can be induced to cause tissue damage while 
keeping the exposures in the overlying tissues below the cavitation threshold. 

Those skilled in the art will appreciate that the invention can be applied for purposes of 
5 tissue ablation, as well as in other procedures where focussed ultrasound is desired. These 
include opening the board-brain barrier, activation of therapeutic agents, occlusion of blood 
vessels, disruption of arteriosclerotic plaques and thrombi, etc. It will also be appreciated that 
the invention can be applied for treatment of humans, rabbits and other animals. The 
embodiments discussed above and shown in the drawings are illustrative only. Other 
10 embodiments, incorporating substitutions, modifications and other changes therein, fall within 
the scope of the invention. These include embodiments with transducer arrays of different sizes, 
shapes and numbers of elements, as well as embodiments with different amplification and 
driving systems. In view of the foregoing, what I claim is: 
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Buchanan et al, 'The Design and Evaluation of an Intracavitary Ultrasound Phased Array for 
Hyperthermia," IEEE Trans. Biomed. Eng. . v. 41, pp. 1178- 1 187 (1994) 
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Design and Experimental Evaluation of 
an Intracavitary Ultrasound Phased 
Array System for Hyperthermia 

Mark T. Buchanan and Kullervo Hynynen 



Abstract— For evaluating the feasibility of treating prostate 
cancer, a 64-element linear ultrasound phased array applicator 
for intracavitary hyperthermia was designed and constructed. 
A 64-channel ultrasound driving system including amplifiers, 
phase shifters, and RF power meters was also developed to drive 
the array. The design of the array and driving equipment are 
presented, as are the results of acoustical held measurements 
and in vino perfused phantom studies performed with Ihe array. 
Several techniques for heating realistically sized tumor volumes 
*ere also investigated. Including single focus scanning and two 
techniques for producing multiple stationary foci. The results 
show that the operation or the array correlated closely with 
the theoretical model. When producing a single stationary focus, 
the array was able to increase tissue temperature by 12°C in 
vitro in perfused phantom. With some minor improvements in 
array design, intracavitary phased arrays could he evaluated in 
a clinical environment. 

I. INTRODUCTION 

INTRACAVITARY ultrasound arrays offer an attractive 
means of inducing local hyperthermia in deep-seated tu- 
mors located near body cavities. By locating the radiators as 
close as possible to the treatment site, problems frequently 
encountered with external Techniques, such as blockage of 
the acoustic window by bone or gas, or simply the inability 
to attain adequate energy penetration, can be avoided. Early 
results using multielement, nonfocuscd arrays of half-cylinder 
transducers operating at 1 .6 MHz suggest that such arrays can 
be clinically useful in the treatment of prostate cancer { I ). 

The proximity of the prostate to the rectum wall makes it 
a good candidate for heating using intracavitary ultrasound 
radiators. Since the prostate is located very near the anus and 
only millimeters away from the rectum wall, an applicator can 
be easily located close to the prostate. The prostate is one of 
the most easily accessible tumor sites, and one that affects a 
large enough population of patients to be potentially clinically 
useful. As such, most of the experiments were conducted with 
the goal of heating the prostate in mind. 

While the nonfocused arrays have shown considerable po- 
tential in healing the prostate, they have two primary Km- 
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nations: the depth to which they can effectively heat, and 
their limited ability to control the power deposition fields 
[2). 131- Both deeper penetration and increased control over 
the power deposition pattern can be achieved using linear 
phased arrays. Since phased arrays can focus their radiated 
energy, they theoretically could heal tissues ro therapeutic 
temperatures deeper than nonfocused arrays, The ultrasound 
power deposition pattern can be electronically tailored as 
necessary to produce therapeutic temperatures within a desired 
volume. In areas whiirc heating would be undesirable, the 
phased array could take advantage of destructive interference 
to minimize power deposition at those locations. The ease 
by which the power dcposiiion pattern can be electronically 
altered in real time provides a means for the compensating for 
varying physiological parameters, patient positioning, and for 
minimizing patient discomfort. 

Several types of ultrasound phased arrays have been 
proposed or built for hyperthermic purposes. These include 
Umemura and Cain's sccior-vonex and concentric ring appli- 
cators |4J, 15]. the cylindrical section applicator developed by 
Ebbini et at. (6). |7J. as well as the tapered array developed 
by Benkescr et al. [8]. Each of these arrays is composed of 
anywhere from 16 to 64 individual elements and operates at 
frequencies between 500 kHz and 150 kHz. While these arrays 
show significant potential, they ate meant to be used in external 
applications and therefore are unsuitable for intracavitary use 
in their reported configurations. 

Previously, a theoretical study on the feasibility of intracav- 
itary phased arrays using half-cylinder elements had been done 
by Diederich and Hynynen [3]. Based on many acoustical and 
thermal simulations, it was concluded thai a practical array 
would be composed of 30 half-cylinder elements with 2.5 
mm center-to-center element spacing operating at 500 kHz. 
The study predicted that the crating lobes formed with this 
array could be negated by using surface cooling. One of 
the major considerations in specifying ihis design was not 
the performance of the array, but the apparent high cost of 
the amplifiers necessary to drive the system. The desire to 
minimize the number of amplifiers led to the specification of 
a 30-element array. 

Initial array designs were based on some recommendations 
from the theoretical study and expanded upon in an attempt 
to construct a practical intracavitary phased array. It became 
apparent that more than 30 amplifiers would be necessary for 
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BUCHANAN AND HYNYNEN: INTRACAVITARY ULTRASOUND PHASED ARRAY SYSTEM 




Fi S . I. Th< 64-c.cmen, arrar and a smgle element showing ihe geomeiry oi .he m * Th, ,™ mounted ro a br„ S shell 19 mm » Mr. ™ 
back or the applicator removable so that vucs could be connected to the inner electrode 



a useful array. This ted to the development of the 64-channcl 
driving system to overcome the limitations imposed by the 
number of amplifiers. This allowed for the evolution of a 6-U 
channel phased array with 1.73 mm centcr-to-cemer spacing 
and a total array lengih of 1 10.5 mm operating at 500 kHz 

II. Materials and Methods 

A. Array Consnuawn 

A 64-element array was constructed using half -cylinder 
,ransducers operating in their resonant radial mode at 500 kHi. 
The array was made by slicing washer- shaped elements with 
a diamond wire saw (Laser Technologies. North Hollywood. 
CA) from 15 mm O.D. by 30-mm long cylinders of PZT- 
4 material (EDO. Salt Lake City. UT^ with silver electrodes 
plated on both the inner and outer wall surfaces. The transducer 
slices were glued together using a silicone adhesive I Dow 
Coming. Midland. MI) with 0.17-mm thick silicon rubber 
spacers^SPC Technology. Chicago. ID between each element 
(producing a dead space of 0.23 mm between elements'). The 
stack of elements was then cut in half along the axis of the 
cylinder and the two half-cylinder sections glued together to 
form the full array. The array was bonded to a brass shell 
to form the complete applicator, as shown in Fig. 1. Wires 
were soldered to the inside wall of each array clement that 
extended the length of the shell to the handle where they were 
each connected to a 2-m RC-178 coaxial cable. Tne electrodes 
on the outer surface of the arTay were electrically connected by 
single wires embedded in beads of silver epoxy (Chomencs. 
Woburn. MA'i running along both edges of the array.. 

5. Driving Hardware 

One of the primary disadvantages of phased anays is the 
increased complexity of the driving equipment. Due to a lock 
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of commercially avaiubie equipment to drive the array. 3 
64-channei ampiiher system was developed. This system is 
compo>ed of phase shifters/duty cycle controllers, amplifiers. 
RF power mcicrs. and impedance matching networks (see 
block diaeram in Fi? 2i and is controlled by an IBM PC 
compatible compuie: a digiUl I/O interface 

Tne amonners arc based on a switching MOSFET design 
ana are designed around inicmauonai Rectifier's IR2M0 dual 
MOSFET onv-r Each of the amplifiers is capable of deliv- 
ering up io 16 W of RF power at 500 kHz into a 50- W load 
from each cnanncl or a soul output power of about 850 W. 
The amplifiers convert digital logic input signals into high 
power sine waves >oiic preserving the phase of the input 
signal The jmpiitude of the output signal is controlled by 
two methods All of me amplifier outputs can be controlled 
simultaneously by adjusting the output voltage on a 1000 W 
DC supply, or the output of each cnannel can be individually 
controlled bv varying the duty cycle of the input signal. The 
result of reducing the duty cycle of the input signal (the duty 
cvcie is the percent of "on" time of the input signal per clock 
■ cycle) is a corresponding decrease in the amplitude of the 
output signal. 

Since the ampiincr: require digital input signals, the phase 
shifting and duiycyeit control is implemented using digital 
ountcrs (101. [HI- These circuits provide 2:.5 3 phase shift 
resolution from 0-360^. and 16 steps of duty cycle (amplitude) 
that allows the output io vary from 0-100% of the output 
power, as allowed by the voltage from the DC supply. While 
the phase shift resolution may seem somewhat limited, studies 
have shown that it is sufficient for the uses described here 

{ 1 The actual RF output power from each channel is monitored 
bv oower meters that measure both the forward and reflected 
RF power 1131. The power meters.. which were also designed 
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/ 


i frequency 


| 500kH2 


c 


1 velocity 


1 SNpm* 


{ a 


j Radius of simple source 


1 V32 


f. 


| wavelength 


! 3 0 mm 



and constructed in house. allow tor verification of (he output 
power on each of the amplifier channels, and allow for 
the monitoring of the output; -o that possibly dangerous 
fault conditions do not occur. B? monitoring the amount of 
reflected power, changes in the electrical matching and the 
anay operation can be easily rri'-mtorrd so that peak array 
performance can be maintained 

Finally, to maximize power tr^ter between the amplifiers 
and array elements, a simple LC :rr. ptdar.ee matching network 
was used on each cnannel bet*t*r. tr.fi oo^er meters and the 
array. This not only maximizes z-. ^z: transfer, but ensures 
greater uniformity m iht accus:;;;! '.ciput signals for a given 
excitation phase and amplitude ir.izc TsperimentallyYv a 
needle hydrophone. Without cie;--:L. marching, the' phase 
shift and amplitude distortion c±.>t: the impedance dif- 
ferences between elements :n Jit ir.i\ vvould have to be 



characterized using a hydrophone and the results used to 
compensate the excitation signals. 

C Actfw.wc Modeling 

The 3-D ultrasound modeling program used *as orieinally 
developed by Diedcrich and Hynynen (3) and is only briefly 
introduced here. The routine models A' cylindrical radiators of 
finite length, radius and separation as each of the ultrasound 
transducers in the anay. and assumes a target media having 
similar acoustical properties as tissue (Tabic I). The routine 
models the surface or each of the cylindrical elements as an 
eveniy spaced grid of simple hemispherical sources and uses 
Huyecns principle to sum the contributions of each source at 
each point in the field. All of the simulations were done with 
the sources spaced on a 1/32 grid. 
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The acoustical pressure field was calculated in the 2-plane The technique solves the Rayleigh-SommerfeJd integral for 
using an array with jV elements, and a field of M control poim> 



Pi*-r) = \ £ 52 ft KOI 



where 



ft(r)a ftp c ^/- 



(2) 



and 7i ? is the number of sources in ihe c -direction. 715 is ihe 
number of sources in the rheta direction. P. is the pressure 
from a sinsie source (Pa). P sa is the pressure amplitude at the 
surface of the source (Pa). r e is the radius of ihe source (m). r 
is the radial distance from the center of the source (m). / is the 
wavelength (m). 6 is the phase of the excitation signal (rad). 
/ is the operating frequency (Hi), and q is the attenuation 
coefficient (Np/ml. 



D. Focusing Techniques 

The single focus case is the simplest form of focusing that 
can be done with a phased array. The single focus is produced 
by setting the phases of the driving signals so that constructive 
interference occurs at the desired focal position. The phase of 
each of the driving signals for an array with ,V elements can be 
calculated from the differences in the path lengths <f. between 
each array element and the focal position by 



4 O 



3) 



where p is the density of the propagation medium, c is the 
speed of sound in the medium, k is the wavenumbcr. S' is the 
surface of the source. u n is the panicle velocity normal to the 
surface of the source, and r m - r n is the distance between «j 
control point ir„) and the surface of a source (r n ). m = I, 

2 >/» and n = 1. 2 A r . 

To simplify programming of this algorithm, (4) can be 
described by 



Hu = p 



(5) 



where u is the excitation vector, p is the complex pressure 
at the control points, and H is the remainder of (4). For 
intracavitary hyperthermia uses, the number of elements in 
the anay (fO is always greater than the number of control 
points iM). This leads to an undcrdetermmed system of 
equations w,th an infinite number of solutions. The minimum 
norm solution lit) can be determined by using a least squares 
approximation of (5) 



(6) 



where /, is the phase (radians) of the 7 th clement. I is the 

wavelength (m). 1 = 1.2 /V and rn = 0. 1,2 The size 

of the focus produced with this technique is usually 100 small 
to heat an entire tumor volume, and therefore other techniques 
must be used to heat larger volumes. 

To heat larger volumes, two techniques were investigated, 
single-focus scanning and mulnple focusing. The single focus 
scanning routine simply stepped a single focus back and forth 
along a predetermined length of the array. The focal depth was 
kept constant while the focus was stepped every 300 ms to the 
next position in the scan. Delivered power was maximal at the 
edges of the scan but was reduced to 649b of maximum power 
at the center of the scan to flatten the temperature distribution 
in the perfused phantom experiments (the power distribution 
was experimentally determined). 

The other technique, multiple focusing, simultaneously pro- 
duces more than one focus within the target volume. The driv- 
ing signals necessary to produce multiple foci were calculated 
using two techniques: split focusing and the pseudo-inverse. 
To create multiple foci with the split focusing technique, the 
array was divided into subarrays. each of which produce a 
single focus in the same manner as previously described. The 
pseudo-inverse method, developed by Ebbini and Cain (15). 
uses a series of control points that represent the magnitude 
of the ultrasound field at given points. A brief summary of 
Ebbini and Cain s technique foUows. 
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where H" * ihe complex conjugate of H. With this technique, 
all of the array elements contribute 10 all of the foci, unlike 
ihe spin focusing technique where each element contributes to 
only one focus. 

111. EXPERIMENTAL TECHNIQUES 

A Acoustical Efficiency Measurements 

The acoustical output of the cylindrical transducer elements 
was measured usmg a radiation force technique (16]. The 
arrav was piaccd in a brass cone with 45 c sides to reflect 
the radialiv emitted ultrasound fields down into the acousucal 
absorber (17). The force on the absorber was measured using 
a Metier AE 160 (Hightstown. NJ) microbalance, while the 
RF power was measured using a Hewlett Packard 438A Rr 
power meter and a Weriatone C2625 (Brewster. KY) dual 
directional coupler. The efficiency was calculated as the ratio 
of the acoustical power to the RF electrical input power. 

B. Ultrasound Fie\d Measurements 

The ultrasound fields were mapped in a tank of degassed, 
deionired water bv mechanically scanning a thermocouple 
embedded in a small f 2-mm diameter) plastic sphere. The ther- 
mocouple was positioned by a three-axis computer controlled 
scannine table." The applicator was mounted on a rotanonai 
device that allowed measurements to be made in a radiai art 
around the arrav bv rotating the array. Measurements were 
made on a 1 x ' 1 mm grid with the recorded data being the 
average of three consecutive measurements. 
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C. In Vitro Kidney Experiments 

Alcohol fixed canine kidneys were used as phantoms for 
studying the heating characteristics of the array. The kidneys 
had previously been prepared as described by Holmes et 
al (18j, ^ were rehydrated prior to use. The experiments 
were conducted at room temperature using degassed, deionized 
water as the perfusate. A metering pump (Fluid-Metering Inc.. 
RHICKC, Oyster Bay, NY) connected to the renal "artery 
circulated water through the kidney while the renal vein was 
allowed to drain into the tank. The kidney was held in place by 
gently sandwiching it between two PVC membranes mounted 
to a Plextglas frame. The applicator was firmly clamped to the 
frame to maintain a fined distance between the surfaces of the 
kidney and array. Fig. 3 shows a diagram of the experimental 
setup. 

Temperatures were measured in the kidney using either 
seven sensor probes sutured in place perpendicular to the 
array through the focal region, or by one or two single sensor 
thermocouples pulled along a path parallel to the array.. In 
experiments conducted with the multiple sensor probes, the 
kidney was exposed to ultrasound for a total of 10 min with 
temperature measurements occurring every 30 s. During each 
temperature measurement, power to the array was disrupted 
for approximately four seconds (one second prior to the 
first reading, and up to three seconds to read all of the 
thermocouples'). Since this technique does not ffive very line 
spatial detail of the temperature distributions, the pull-back 
technique was more frequently used. 

The pull-back experiments were conducted by pulling one 
or two single uncoated thermocouples (0.05-mm wirei by \\ 
computer controlled stepper motor along a track parallel 10 
the array through the kidney tissue. The Temperatures were 
measured every 1 mm and the average of three readings was 
recorded. Prior to each of the experiments, a baseline tem- 
perature was established along the path of the thermocouples. 
The kidney was exposed to ultrasound for 20 min 10 allow 
the temperatures to reach steady-state before the temperature 
profiles were measured. The difference between the baseline 
and final measurement was used to calculate the temperature 
rise. The kidney was allowed to cool to room temperature 
(typically, 30 min) before to the next experiment was started. 

The thermocouples were mainly located in the medulla of 
the kidney. The steady-state temperatures in the medulla have 
been shown to be a strong function of the flow into the kidney 
[19). The flow values were kept relatively low in order to 
simulate the perfusion in the prostate (20). [21]. 

IV. Results 

A. Acoustical Efficiency Measurements 

The acoustical efficiency of the 500-kHz PZT-4 transducer 
material was measured for a full cylinder. 30 mm Ions, and 
for two half-cylinder arrays wiih 2.5 mm and l.S-mm center- 
to-center element spacing (2.2 mm and 1.5 mm element 
widths) with total array lengths of 24." mm and 2S.5 mm. 
respectively. Small arrays were tested so that the resulting 
field would be collimatcd. unlike a single array eiemeni 




Fig. ?. Diaer»n of the »»> xim, kidney experimental »etup 
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The results of these experiments are sho^n in Fie 4. None 
of the arrays exhibited very high efficiencies despite the low 
operating frequency Though not shown here, some results 
were verified using ca'iohmeinc techniques. While the 30-mm 
long full cylinder element had an efficiency of 71 %. the arrays 
exhibited efficiencies of 27% and )J% for the 2.2-mm and 1.5- 
mm arrays, respectively. These efficiencies indicate mat about 
809c of the electrical power delivered to the ).5-mrn wide 
elements used in the phased arrays is lost either a? heat in the 
transducer or in the eiectncal matching and transmission lines 
between the amplifier and the array. 

B. Ultrasound Field Measurements 

The effects of element spacing were studied by measuring 
the ultrasound field distributions produced by two 16-eiemeni 
arrays with l.^mm and 2.5-mm center-ro-cemcr element 
spacing Fig. 5t a > and tb) shows the acousticai field plots made 
in a plane paraiiel to the an-ay and normal to the surface of 
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the hail cylinder elements for two 16-element arrays with t.S- 
mm and 2.5-mm element spacing along with their respective 
simulated fields (Fig. 5(c) and (d)). Both array? were focused 
30 mm from the array surface and 15 mm above the central 
axis of the array. Tne most significant effect illustrated here is 
the generation of the grating lobe by the 16-element array with 
2.5-mm center-to-center element spacing (Fig. 5(b) and (d)). 
Since such a large grating lobe not only reduces the power 
of the focus but also could cause heating in unwanted areas, 
arrays must be designed to minimize grating lobe formation. 
The array with l.S-mm cemer-to-cemer spacing produced 
much better field distributions than the array with 2.5-mm 
cemer-to-cemer spacing. Therefore, in order to have center- 
lo-cemer spacing 1.8 mm or smaller. 1.5-mm wide elements 
were used in the final array design and the dead space between 
elements was reduced to 0.23 mm. 

The complete measured and simulated plots of the ultra- 
sound field produced by the 64-element array with 1.73-mm 
element spacing axe shown in Fig. 6(a) and (b). respectively. 
The array was focused 40 mm deep, and 30 mm from the 
central axis. Note thai only a small grating lobe is produced 
by this array even though it is focused a substantial distance 
from its central axis. In fact, this array is capable of focusing to 
about 35 c from the central axis of the array without producing 
significant grating lobes. 

29 



So far. all of the acoustical field measurements presented 
have been tn a plane parallel to the length of the array and 
normal to the surface of the array. Ideally, the acousticai field 
would taper on toward the edges of the array but be more or 
less uniform around the arc of the an-ay. Unfortunately, this 
ideal docs not match reality, as can be seen in Fig. 7. The 
acoustic field was measured as a function o: rotation angle 
along a tine by fixing the position of the thermocouple in 
the focal region and rotating the array about its axis. The 
acoustical imensiry is not at all uniform and. in fact, vanes 
as much as 50% before tapering off at the edges. AH of the 
2-D acoustical field plots shown here were made at the- 0° 
rotation angle where the intensity is only about 50% of the 
peak. Measurements made, but not shown here, show that this 
fluctuation in pressure amplitude due to the rotation angle only 
alters the peak intensity and does not effect the overall shape 
of the acoustic field. 

Fig. Si a) and (b) shows the acoustical field plois produced 
usine the cwo multiple focus techniques: pseudo-inverse and 
split focusing. Below them, in Fig. S(c) and (d). are their 
respective simulation results. Both techniques were used to 
synthesize foci 40 mm from the array surface and 10 mm 
on either side of the central axis of the array (20 mm apart) 
using the 64-element array. Notice that the pseudo- inverse 
technique produces sharper foci than those produced by the 
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splii focusing method. This occurs because ihe pseudo- inverse 
technique utilizes the entire length of the array allowing it to 
produce sharper foci than the split focusing technique which, 
since it divides the array into two subarrays, effectively uses 
array lengths only half as long as the actual length of the array. 

C. In Vitro Kidney Experiments 

Fig. 9(a) shows the temperature nse ver^u* time along a 
fixed seven-sensor thermocouple probe located perpendicular 
to the array. The distances marked denote the distance of 
the thermocouple from the edge of the kidney nearest to 
the array. The array, with 32 active elements, was focused 
along its centrai axis and 10 mm into the kidney. (Thirty- 
two active elements were used since the array length in this 
configuration could adequately cover the kidney.) The kidney 
was perfused at a rate of 2.9 kg m _J s~ l . Note that the kidney 
is only 35 mm thick and thus, the temperatures close to the 
surfaces of the kidney are dominated by cooling caused by 



the circulating water in the water bath. All of the following 
results were obtained in the center of the kidney where the 
focus was located and the surface effects were smallest. The 
temperature rises achieved in the middle of the kidney indicare 
thar therapeutic temperatures can be achieved with a stationary 
single focus at realistic perfusion levels. Typical cross-section 
profile^ measured in the middle of the kidney with the pull- 
back technique tor a center focus and 15 mm on either side 
are shown in Fig. 9(b). 

The three different techniques for heating larger volumes: 
pseudo-inverse, split focusing, and scanned single foci arc 
compared in Fig. I Of a ) and (b). The two multiple focusing 
techniques were used to produce foci 40 mm from the surface 
of the array and 10 mm on either side of the array s central 
aMS, while the scanning routine produced a series of single 
foci 40 mm from the array surface with a total span of 
20 mm. The thermocouples were located 8 mm and 15 
mm from the surface of the kidney nearest to the array 
(35 mm and 42 mm from the surface of the array). The 
scanning 'technique produced a narrower profile than the two 
multiple focusing techniques and the distribution lacks the 
lemperature drop between foci. The two multiple focusing 
techniques produce virtually identical profiles, although the 
pseudo- inverse technique does not produce as large of a 
lemperature rise on the deeper thermocouple (beyond the 
focus > as does the split focusing technique. This IS caused b> 
the increased sharpness of the foci produced by the pseudo- 
inverse technique, as previously explained. 



v. Discussion and Summary" 
An intracavitary ultrasound phased array composed of half- 
cylinder transducer elements has been constructed for inducing 
hyperthermia in the prostate via the rectum. A 64- channel 
amplifier system has also been designed and constructed to 
drive the phased array. As was Shown by the acoustical field 
plots, the array ts capable of producing focused fields as well as 
fields containing more than one focus. The array is current!) 
capable of producing a 12°C temperature nse in a perfused 
phantom using a stationary focus, and smaller temperature 
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increases were achieved in larger volumes by scanning the 
focus or by creating multiple stationary foci. This was achieved 
with flows to the kidney that should simulate the relatively low 
perfusion in prostate [20), [21]. 

What has yet to be conclusively shown is such an array's 
ability to heat tissues to therapeutic tempera lures in a regu- 
lar and predictable manner. The peak temperatures achieved 
within the perfused phantoms tended to vary considerably from 
experiment to experiment, making direct comparisons between 
the absolute temperature profiles difficult. Therefore, the tem- 
perature profiles that compare various focusing techniques 
were normalized to better illustrate the overall differences 
in the shape of the temperature profiles. The fluctuations in 
temperature were caused by a variety of problems, including 
morphological differences in kidneys and the location of the 
thermocouples within the kidneys, the inefficiency of the array, 
and the open loop manner in which the array was coerated. 

The electrical efficiency of the array would usually drop 
considerably during the first experiment, and during each sub- 
sequent experiment due to changes in the electrical impedance 
of the array elements. Two primary factors weTe responsible 
for the observed changes in the electrical impedance of the 
array elements: The first was a thermally induced impedance 
drift caused by the array self-heating during sonication. The 
second impedance shift that occurred from experiment to 



experiment because the impedance of the array elements did 
not return to their original impedance after cooling. The 
impedance shifts made it difficult to accurately control the 
radiated power delivered to the kidney. Other problems with 
the array included arcing between elements and water seepage 
behind the array that reduced the array efficiency in subsequent 
experiments until the array was repaired and remaiched. With 
a more carefully constructed array, most of the power limiting 
problems experienced can be avoided. 

The acoustical field plots were in good agreement with the 
theory, and the techniques for heating larger volumes were 
all functional, though the multiple focusing techniques were 
the most effective. Overall, phased arrays show considerable 
potential for improvement over currently used intracavitary 
ultrasound hyperthermia system. 

While phased arrays allow significantly more control over 
the acoustical field, the current design using half-cylinder 
radiators still lacks control in the angular direction (around the 
arc of the array). Additionally, since the cylindrical radiators 
do not have uniform angular intensities, the angular heating 
pattern is somewhat degraded, though thermal conduction will 
probably smooth the resulting temperature distribution. Similar 
fluctuations in the angular field disrributions have been shown 
for other cylindrical transducers (2) T [11). Radial control 
could be achieved by dividing the half-cylinder elements into 
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pie-shaped subelements and driving each subelement indepen- 
dently. However, this would not only dramatically increase 
the number of amplifiers, but the resulting subelements would 
have virtually no inner electrode due to the wall thickness 
of the element. This would seriously degrade the already low 
acoustical efficiency of the cylindrical radiators. 

A possibly better array design would utilize planer array 
elements mounted on a routing platform. By tilting the array 
back and forth, the same volume could be treated as with 
the cylindrical radiators. By controlling the acoustic field and 
output power as a function of tilt angle. 2-D control over 
the power deposition held can be achieved without requiring 
additional ampliners. An additional improvement on the array 
design would be to increase the operating frequency. While the 
500-kHz operating frequency allowed larger element sizes and 



Fig 10 Nornuiiico tcmperature/po»*er profiles measureo m a perfused 
kidnev iff mumpie t'ocuMng u«mg ine pseudo-inverse and scSu focusing 
technique*. »» a* smjlc-focus ^canning The multiple loci *ire created 
40 mm oeep. 10 mm on either side ol the central axis while me seanntng 
vj* done «.nh « >mek focu< 40 mm deep, and scanned 10 mm on either 
vide of the ecnirti an> Tne measurements «<crc made 6 mm (a) and 15 mm 
(bl from the surface or" ihe kidney with the array positioned 27 mm from the 
kidnev surface 



minimized grating lobe formation, higher frequencies would 
increase the po*er absorption in tissues. Since most tumor 
sites tor which this applicator would be used occur near 
the cavity v. all. the deep penetration allowed by the 500- 
kHz operating frequency is not necessary. The difficulty in 
increasing the operating frequency is that the element size and 
center-io-cemer eicmem spacing would need to decrease as 
the frequency increased, making array design more difficult. 

As a conclusion, the intracavitary, electrically focused array 
did demonstrate, at practical depths for prostate treatments 
(about 30 mm), the feasibility of using phased arrays for 
intracavitary hyperthermia purposes. With a more carefully 
constructed aray. most of the power limiting problems ex- 
perienced can be avoided. The acoustical field plots were 
in good agreement wiih the theory, and the techniques for 
heating larger volumes were all functional, though the multiple 
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focusing techniques were the most effective. Overall phased 
arrays show considerable improvement over the currently used 
intracavitary ultrasound hyperthermia system [1). However, it 
does not offer power control in the angular direction (around 
the axis of the applicator) and thus, new designs allowing 3-D 
control of the power field should be developed before clinical 
testing. 
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Focused high-power ultrasound beams are well suited for noninvasive local destruction of deep 
target volumes. In order to avoid cavitation and to utilize only thermal tissue damage. ,high fre- 
quencies (1-5 MHi) are used in ultrasonic surgery. However, the focal spots generated by sharply 
focused transducers become so small that only small tumors can be treated in a reasonable time. 
Phased array ultrasound transducers can be employed to electronically scan a focal spot or to 
produce multiple foci in the desired region to increase the treated volume. In this article, theoretical 
and experimental studies of spherically curved square-element phased arrays for use in ultrasonic 
surgery wctc performed. The simulation results were compared with experimental results from a 
16-eiement array. It was shown that the phased array could control the necrosed tissue volume by 
using closely spaced multiple foci. The phased array can also be used to enlarge a necrosed tissue 
volume in only one direction at a time. i.e.. lateral or longitudinal The spherically curved 16 
square -element phased array can produce useful results by varying the phase and amplitude setting. 
Four focal points can be easily generated with a distance of two or four wavelengths between the 
two closest peaks. The maximum necrosed tissue volume generated by the array can be up to 
sixteen times the volume induced by a similar spherical transducer. Therefore the treatment time 
could be reduced compared with single transducer treatment. 



Key words: phased array, ultrasonic surgery 
I INTRODUCTION 

In ultrasound surgery, to avoid cavitation and have a sharp 
boundary between the rumor and normal tissue, high fre- 
quency focused transducers have to be used. 1 The focal spots 
generated by sharply focused ultrasound transducers are 
small when compared with the diameter of many tumors. 
Obviously this type of focus is not efficient to treat large 
tumors, which would require a large number of exposures to 
cover the whole target volume. However, in order to accu- 
racy treat the target volume close to critical structures, 
s . focal spots may be required during pan of the therapy. 
Thus, controllable focal spot si2e is required. The most at- 
tractive method to obtain control over the focal spot sue is to 
use phased arrays. 

Phased array applicators were introduced Id ultrasound 
hyperthermia cancer therapy in the early 19S0's. During the 
past decade, many efforts have been made to investigate the 
advantages of phased arrays in hyperthermia, and several 
phased array applicators have been developed. Phased array 
applicators can be divided into the following categories: an- 
nular or concentne-ring arrays, 2,3 stacked linear-phased 
arrays, 4 sector-vonex arrays^ tapered linear-phased .irrays. 5 
cylindrical-section arrays, 6 and square-element spherical- 
section arrays. ' It was shown that ultrasound phased arrays 
can provide good control over the heating pattern with flex- 
ibility in moving the focus and producing multiple foci. Pre- 
vious research in ultrasound phased arrays has been concen- 
trated in hypenhermia cancer therapy. Relatively low 
frequencies (about 0.5 MHz) and small element size com- 
pared to the wavelength) were employed in these studies 
The studies on square-element spherical- section arrays em- 
phasized using a large number of elements, up to a few hun- 



dred. The main disadvantage associated with using a large 
number of elements is that the same number of amplifiers 
and electronic circuits are also required. 

The motivation of this article was to demonstrate that a 
phased array system with a small number of elements can be 
utilized 10 provide control over the focal spot size for surgi- 
cal purposes. Tne excitation signals of the elements vary in 
both amplitude and phase. Hence the major task in utilizing 
phased arrays is to determine the phases and amplitudes 
needed 10 produce desired ultrasound fields. The inverse 
technique introduced by Ebbini and Cain 6 can be used for 
direct synthesis of ultrasound fields with multiple foci. Sev- 
eral amplitude and phase settings were calculated for differ- 
ent sets of selected control points. The limitation of the num- 
ber of utilizable control points was also investigated for a 
given array. Several simplified amplitude and phase settings 
based on the calculated amplitudes and phases were em- 
ployed for ultrasound field calculations. These driving signal 
sets can be utilized, when different focal spot sizes arc re- 
quired for the array proposed here. The transient bioheat 
transfer equation was employed to estimate die temperature 
elevation due to the ultrasound power deposition. Then the 
necrosed tissue volume was predicted by the isothermal dose 
volume. Computer programs were aiso used to do a paramet- 
ric study to invesueate the influence of the dimensions and 
frequency of the array on the necrosed tissue volume. 

II. MATERIALS AND METHODS 

A. Square-element array with spherical surface 

The phased array was made from a spherically curved 
transducer which was divided into small square elements. 
The whole array was airbacked. It was constructed from a 
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FiG. I. Configuration of *.r.t «quare-cremcnt sphenc&llv curved pr.:t.<eO way 
with i« coordinate ?>»icn>. 



where j ~ n/-T. p is the density of medium, c is the ipeed 
of sound, k is the real wave number, o n = A n e ,tvt * #J is the 
complex excitation source of the nth element with amplitude 
A r and phase B n . r m , is the distance frorr. a point (x fl .>** ,z r ) 
on the nth element to the field point of interest {x„ >y m .z m ), 
5„ is the area of the nth element. 

In ultrasonic surgery, waicr is used as a coupling medium 
between the transduce?* and tissue interface. Water is consid- 
ered a nonancnuatmg^mediunx and tissue a weakly attenuat- 
ing medium. The effects of anenuauon can be satisfactorily 
described by replacing c~ Ji '"<* wuh e m expression 
(1 1, where k r = k-ja is the complex wave number with at- 
tenuanon coefficient a." 12 

The power deposition for the desired volume :s given by 

|p«u.y.ci| 3 



Spherically curved PZT4 bowl, having a diameter of 100 
mm. a radius of curvature of 130 mm. and a frequency of 1 .4 
MHz. The bowl was cut into 16 elements, each with a length 
of 20 mm per side. A 0.3-mm space between the elements 
was filled with silicone rubber for electrical and mechanical 
isolation. Each or the elements was connected to an LC 
matching circuit to match the impedance to 50 O. and 0 e . The 
array was driven *ith a custom made 16 channel amplifier 
(Labthermics. Champaign. Illinois). The phase and ampli- 
tude were controlled by RF signals feeding the amplifiers 
The driving signals ^ere generated by an ui-house manufac- 
tured digital circuit.' The amplitude and phase of each input 
signal were digital i> controlled with 3 -bit resolution in am- 
plitude and 4-bu resolution in phase. The number of ck- 
mems in each row ano each column were the same so thai 
the whole applicator u as a square shaped focused Transducer. 
Tnc configuration 31 the phased array for the experiments •> 
shown in Fig. I . 

B. Ultrasound field measurements 

The relative pressure amplitude squared distributions 
were measured in degassed waier using a needle hydrophone 
(active spot size 1 mmi scanned across the focal region. The 
needle hydrophone was moved by siepper motors, typically 
with 0.1 -mm steps across the beam. The total acoustic power 
was measured ujm: a radiation force technique. 

C. Computation of ultrasonic fields 

Consider an ultrasound phased amy with N eiemenu. 
employed to produce an ultrasonic field in a nonattenuaung 
medium. Assume me coordinate system is defined as in Fig. 
I. According to the theory of ultrasonic radiators developed 
by O'Neil" 1 and the principle of superposition of acoustic 
pressure, the pre«ure field due to this transducer can be 
evaluated by ihe'Rasleirh-Sommerfeld integral The inte- 
gral is 



jpek v 



12) 



where M»bi ,s ,nc absorption coefficient and Z-pc is the 
impedance of the medium. 



D. Inverse tochniqua 

The inverse lechnique can be used to calculate the ampli- 
tude anc ona^e senmgs from selectee control points wnere 
the oestrec pressure values are given. Defining 

"•' ,.; p r w „ 

tnen lor M field points, expression \ 1 f can cc written :n 
matnx notation as 

P-HL u ' 

wnere 



; P\ 

,! P: 
i 

\ 



H- 



H 

i 



The elements of matrix H are evaluated by numerical inte- 
gration usinc expression (3). Equation (41 has two important 
features. First, if the H matrix is calculated ar.c saved for a 
desired neid. then the pressure field can be evaluated using 
Eq. i4 1 instead of Eq. (1) for a given excitation source marrix 
U. Computation time can be saved by using :his method 
when varying the excitation source over me same calculation 
volume. Second, tor a oestred pressure field p2t:em. i.e.. tor 
a given matrix P. '.he excitation source U car. ce calculated 
h> 

U-H'P 

wnere H" it the pseudoinwrse matrix of H. Usually the tout 
number of control points arc much less than the total number 
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Table I. Selected comrol poirus u«d for im-ewe takul^.ous for 16 *quare- 
■' dement spherical curved pha>ed amy. The ccmirol points are locked ai 
* :« 129 nail plane. Note ih« units for .« and r are millimeter*. 

v Ph«*e 



of «cii2Uon sources in this case. The W matrix can be 
foi-*d by utilizing single value decomposition. 

.c control points were selected in a plane ai a desired 
focal position. All the points were evenly distributed on a 
circle. Since the pressure at a control point is a complex 
number, both the magnitude and the phase are required to 
perform the inverse calculations. In selecting control points 
in this manner, it is natural to make the magnitudes the same 
for all the points, and the phases either in phase or out of 
phase so that the phase rotates around the central axis. The 
latter setup gives destructive interference on the central axis 
thus, eliminating potential hot spots on the axis/ The se- 
lected control points used in the inverse calculations ar; 
given in Table 1. 

E. Thermal modeling 

An approximate temperature response to the power depo- 
sition was predicted by the transient bioheat transfer equa- 
tion. The differential equation is 

15) 



where T is the temperature Hat Emei . . H_$£&X3iti» HX^iL 
p, is the density of the tissue, c, is the specific heat of the 
tissue. k, is the ihermal conductivity of the tissue, h- is the 
blood perfusion rat, . c h is th. specific heal of the bloc J. 7"„ 
is the arterial blood temperature, and £(.t.v.-) w the acous- 
tic power deposition rate per unit mass. The ihermal response 
was simulated in a homogeneous medium. A surface tem- 
perature of 37 C C on the cube and an initial temperature of 
37 *C inside the cube were used as boundary and initial con- 
ditions for all of the computations.^ A numerical finite differ- 
ence method was employed to solvit Eq. (5). Previous studies 
have shown :hat there are sever*! parameters which affect 
the temperature elevation. 1 * Temperature elevation has been 
shown to be almost independent of (he blood perfusion rate 
for short ultrasound pulses. 14 For this reason, ultrasound 
pulse durations of 1, 5. ard 10 s were used in this study. The 
maximum temperature reached in all simulations was kept 
the same (SO e O by adjusting the input power 

F. Thermal c*ose calculation end necrosed tissue 
volume estimation 

The thermal dose calculation was based on the technique 
suggested by Sapareto and Dewey. 15 Using this technique, 
the accumulated thermal dose was calculated at a reference 
temperature by numencal integration under different tem- 
perature profiles. The thermal dose. i.e.. equivalent time, at 
the reference temperature can be evaluated by 

DoseUr ()= p" rff-"£V-- r *Ai. 

(61 

where T rr , is [he reference temperature, f*^," 'fteitws+'roohnj 
is the final time. A; is a small time interval. T 2 , is the aver- 
ase icmperaiure dunnc time Af. and R is a parameter giver. 

by 

'0.5. if fms*43 C C 

* = \ 0.25. oiherwise 
The necrosed tissue volume is estimated by the isothermal 
dose volume surrounded for 240 nun at the reference tem- 
perature of 43 °C This technique has been found to be 2 
reasonable model tor predicting tissue necroses induced by a 
spherically curbed transducer. 1 * 

A two iavered medium water-tissue was assumed in the 
simulanons. The speed of sound and the density were 1500 
rrvs and 998 kg/nr. respectively, for both media. The attenu- 
ation coefficient of the tissue was assumed to be 10 Np/nV 
MHz. The thermal properties of the tissue are given in Tabic 
II. 



HI. RESULTS 

A. Comparison of experimental results with 
simulations 

It was necessarv to verify the numencal model using the 
16 square-element' phased array with a spherical surface be- 
fore relying on the simulations. The experimental and simu- 
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kpm 1 
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lated profiles of the pressure amplirude squared at the acous- 
tic focus are displayed in Fig. 2. AH the data were 
normalized to 1 by dividing by the maximum pressure am- 
plirude squared for each curve. Figure 2(a) shows the results 
for the case with uniform excitation sources. There was good 
agreement between the simulations and experimental results 
for the main beam. When the array was excited by various 
phase settings, the simulations reasonably predicted the side- 
lobe distributions found in the experiments [Figs, 2(b)— 2(d)] 
The maximum power output measured was at least 14 W for 
each transducer element. 

B. Simulation results 

The phased array ultrasound transducer used in the ex- 
periments was simulated first. A calculation volume of 30 
X 30X90 mm 3 was used for this phased array. The axial dis- 
tance from the center of the transducer to the water-tissue 
interface was 70 mm. Tne inverse technique results are 
shown tirst. The amplitude and phase settings calculated by 
the inverse technique corresponding to selected control 
points arc given in Tabic III (only the first four are shown. 
Tne contour plots o: power deposition calculated for these 
array settings are displayed in Fig. 3. For tour control points 
only one focus was produced, which was centered on tne 
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axis when the pressure at all of the control points were in 
phase [Figs. 3 rial and 3(b)]. Four focal spots were obtained if 
the pressure at the control points had 90.0" phase shifts from 
one point to another [Figs. 3fc) and 3(d)] Comparing Figs. 
3ic> and 3tdi with Figs. 3iei and 3(0, the pea* to peak dis- 
tance was increased when the control points on a 1.25-mm 
ramus circle increased to a 3.25-mm radius circle. However, 
if tne cor.uoi points were on z 5.25-mm radius cucle. the 
local spots were at the same location as Figs. 3(ci and 3(d). 
Tne nejd distributions mot shown here* were very similar to 
Figs 3 cj and 3«'dl For six control points. s:x focal spots 
were obtainec at control locations with different partcms 
[Figs. 3 i and 3(h) j. ana two additional focal spots were 2t 
the center. When the number of control points increased to 
eight, the power deposition pattern had the maximum pres- 
sure amplitude ;n front of the focal plane with strong peaks 
at the water- tissue intenace [Figs. 3(i» and 3u<l These re- 
sults were used to select simplified amplitude and phase set- 
tings wnich *cre then used in necrosed tissue voiume simu- 
lations (Table IV). 

The isothermal doses for tissue necrosis for five amplitude 
and phase settings are Displayed in Fig. 4. Let us consider the 
ultrasound pulse* duration 10-s case (solid line-. For the uni- 
form excitation case, the focus was located at a depth of 59 
mm from the intenace. Based on the 240-min isothermal 
dose [Figs. 4 (a . and 4(b)]. the shape of the necroseo tissue 
voiume was ciose to an ellipsoid. Tne calculated volume was 
about S5 mm : \ measuring 3 mm laterally and 18 mm longi- 
tudinally. Isothermal doses are also shown in the same ns- 
ures for the ultrasound pulse durations of 5 and 1 s. where 
the input power was admstcd so that the maximum tempera- 
ture was kept at SO c C>or case 1 of Table IV. the necrosed 
tissue length [Fig 4. cm was 42 mm. and the width [Fig. 
4(di] was 3.4 mm. The length was more man double and the 
width *as slightly enlarged compared to the uniform excita- 
tion case. The computed" nec-os'C* tissue volume was about 
339 mm\ Ficures 4ici and 4(f) snow case II of Table IV. for 
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Ta8L£ IV, Tne simplified phase and amplitude settings for the 16 sqvar*- 
• element spherically cwved phased amy uses' m ihe necrosed liuue volume 
calculations. 
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which the necrosed tissue length and width were 20.8 and 
5.1 mm, respectively. The length was slightly longer and the 
width was almost double compared with the uniform excita- 
tion case. Notice that this case is equivalent to the situation 
in which the phased array had only four elements with sides 
of 39.4 mm. The calculated necrosed tissue volume was 
about 377 mm ; . Figures 4(g) and 4(h) indicate that there 
were four focai points in the power deposition for case III of 
Table IV. The shortest peak to peak distance was about 4V 
where X is the wavelength. The whole volume was the sum- 
mation of the four individual smali necrosed tissue volumes 
which were produced by the four focai points. Each small 
volume was about the same size as the volume in the uni- 
form excitation case for the 1- and 5-s sonicauons. For the 
10-s sonication a united volume was created wnh the total 
voiume of 810 mm*. For the amplitude and phase setuns of 
case IV in Table IV. the power deposition (not shown here? 
indicated that mere were four strong focai points \ depth 54 
mm) surrounded by four small focal points. The shortest dis- 
tance between the strong peaW was about 2k. Simultaneous 
focusing at these points was performed to enlarge tne heated 
volume compared to the uniform excitation case. The calcu- 
lated isothermal dose [Figs. 4u) and 4'j)] from the tempera- 
ture distribution had a maximum length of ?S.5 mm and i 
maximum width of 7.2 mm. and the necrosed tissue volume 
was close to 1 .4 cm 5 . 

To illustrate the effect of frequency on necrosed tissue 
voiume production, the isothermal doses for case IV of Table 
IV with an ultrasound puise duration 10 s axe shown in Fi2S. 
5(a) and 5(b) The calculated volume was 60X60x90 mnv 
for 0.5 MHz. 40X40X90 mm* tor 1.0 MHz. and 30x30x90 
mm* for 1.4 and 2.0 MHz. The distance from the transducer 
to tne interface was kept the same {70 mm'i. while the frt- 
quency wax varied. The ratio of the necrosed tissue length to 
widtn was almost independent of frequency. P.ie size in- 
creased as the rreque^y decreased. Tne shapes were simiiar 
for freouencies of 1 ic 2 MHz. 
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To understand the effect of radius of curvature on the 
necrosed tissue volume, the isothermal doses for case IV of 
Table IV. with an ultrasound pulse'duration of 10 s and op- 
crating frequency of 1.4 MHz are displayed in Figs. 5(c) and 
5(d). The distance from the transducer to the interface was 
varied according to the radius of curvature to maintain a 
constant focal depth. Tne necrosed tissue volume increased 
rapidly along the longitudinal direction toward the interface 
as the radius of curvature increased. However, the diameter 
of the necrosed ussue volume increased only slightly in the 
focal plane. 

To show me capability of the phased a/ray to move the 
tocus. contour plots of the power deposition for tne maxi- 
mum displacement arc displayed in Fig. 6. Tne simulation 
model *as tne same as the one used to generate F;g. 3. The 
phase distributions were obtained based on direct calcula- 
tions. : c.. phase e^lirdJK, n = 1 ,2 ,V. where d> is the 

distance from ihe center of the element to a desired focus. 
Figure 6iaj snows the results for the case with uniform am- 
pliruoe and pnase setting. The focus could be moved 7 mm 
closer to me interface [Fig. 6(b)). or 7 mm deeper [Fig. 6(c)] 
compared w,th Fig. 6'a). Figure 6(d) shows the phased array 
focus shitted 1 5 mm off the central axis. When iareer dis- 
placements a ere attempted, the phase increment between ad- 
jacent elements exceeded tt/2. Distributions generated by 
those Dnascs are no longer single, strongly focused ultra- 
sound r.z.z: Therefore these were the maximum displace 
menis arnievabie with this array in ihe sense of keeping the 
fields :n ihe region of focus similar to those for the uniform 
excitation case 

Finally. r\:. 7 sho^s the input power requirements for the 
uniform c\c;:at:on and four amplitude and phase semnss 
gner. ir. 7a::. c iV For uiirasound pulse durations cf 5 or 10 
S. tne po* e: ranged from 30 to 250 W For a i -s pulse, the 
input pouer could be as high as 900 W. depenaing on the 
ampin-jar sr.c phaie settings. 

IV DISCUSSION 

The experimental and simulation results showed that a 
I6-eierr.cn: pnased array can offer significant control o%cr 
the size and shape of ihe necrosed tissue voiume during u!- 
trasouno surgery. The experiments showed that sucn an arra% 
can be ccnsiruc.ed and it can deliver enough power for sur- 
gical purposes. Tne simulated and experimental ultrasound 
field distnoutions agreed reasonably well. The differences 
may De cuz :c uneven eiement size and power output in the 
phased arrav as wcii as measurement en*ors. Since the peak 
to peak disiance was only a few miUimeters. ihe ultrasound 
detector may contribute some measurement errors. The simu- 
lation mode; accurately predicted the locations of sidelobes 
and mam beams except that it over predicted the magnitude 
of the sidelobes When adjacent elements of a phased array 
are verv ciose to each otner. the mutual coupling between the 
elements would have an effect on the ultrasound f.eld. The 
difference oetween tne sidelobes in the experiments and the 
simulations may be due to the lack of such crosstalk in the 
calculations 

The auove results proved that the large eiement spi^n- 
cally curved phased array can enlarge the necrosed tissue 
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FiO. 5. The isothermal dose for cue IV of Table IV unh various frequencies (a), ibl. and *«tn *anous ndu of curM-.ure tc). fdl The ultrasourui duI* dur*xion 
used in the themul dose calculation* wis 10 s. The distance inm the transducer to me menace i 50. iOC. 7S and 22 r.m icr raou of curvature of 20C 
160. 130. and SO mm. respectively (ah Icl is axial plane distnouticn. and ib>. idi is. me r'ccai p;*n; amnbunor. 



volume. The phased array can also enlarge ihe necrosed tis- 
sue volume in only one direction at a time, if desired. It is 
important to be able to control the focal spot size so thai 2 
large tumor could be treated in a reasonable time. The con- 
struction of the whole system was relatively simple due 10 
the small number of elements. The 16-elemcnt phased array 
can aiso shift the focus off the central axis or move the focus 
along the central axis. However the distance is limited to 
= 1.5 mm laterally and 14 mm in the axial direction. In mov- 
ing the focus along the central axis, the array is similar to 2 
concentric-ring anay with two rings. Theoretically, the maxi- 
mum phase increment between adjacent elements is rr 
Therefore, the maximum possible phase difference between 
the smallest and largesr phases is tt when moving the focus 
along the central axis. From geometric considerations, this 
phase difference produces displacements along ihe central 
axis of up to 23 mm C10 mm closer, 13 mm deeper). In 
shifting the focus sideways, it is similar to a cylindrical- 
secuon array with four elements. The maximum possible 
phase difference between the smallest and largest phases is 
3* for shining.the focus sideways. Geometrically, this phase 
difference shifts the focus =3.4 mm laterally. However, the 
simulations snowed that the phase increment between the 
adjacent elements shot'' 1 be less th*n tt/2 10 generate a 
single strongly focused ultrasound field. 
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The 16-eiemen: phased array can generate icur focal 
points with a pea* to peak distance as short as two wave- 
lengths The maximum distance between the closest peaks is 
limited to arout four wavelengths. When the four control 
points were on z circle of radius 5.25 mm. the phase incre- 
ment between idiacent elements exceeded tt for this array 
Therefore. 11 :s physically impossible to use the phase and 
amplitude settine obtained by inverse techniques to produce 
the four tocai spots at :hs control points. Although Eq (4) 
represents an unoerdetermind system, it produced a solution 
at the control Domts When multiple foci are separated by 
distances larger than 4.6 mm, the necrosed tissue volume 
becomes a tew individual smai':er volumes, which arc not 
united for shon ultrasound puise durauons 

Bv decreasing the frequency, the necrosed tissue volune 
can be eniarced'because the focal spoi increases due to the 
increased wavelength. But the ratio of the necrosed tissue 
length to the width was kept almost the same. As the radius 
of curvature is increased, the necrosed tissue volume can also 
be enlarged. The necrosed tissue volume is increasefl mainly 
m the axial direction. This agrees with previous experience 
using single clement spherically curved transducers. * 

Phascd arrays require more input power than similar 
vncle focused transmicrrs >o reacn the same temperature 
level for the same ultrasonic puhe duration due 10 increased 
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focal spot size. The tesi array measurements showed that 
thcrequired power levels can be generated in practice at least 
V and 10-s exposures and thus, the proposed technique is 
fusible. 

In Order to obtain the desired field partem, several field 
control points are necessary to perform the inverse calcula- 
tions. However, the number of uiiliiable field control points 
is limited by the number of phased array elements. Hence the 
number of field patterns that produce significantly different 
shapes or sizes of the measured volume is limited. Although 
only few amplitude and phase settings were presented in this 
study, the results provided information which could aid in 
planning ultrasonic surgery using a phased array system with 
a small number of elements. These results have not been 
optimized: nevertheless they illustrate the feasibility and 
range of focal size and shape obtainable with a 16-element 
array. The calculations and experiments can be used to obtain 
a set of different focal spots for a given array. This limited 
number of focal spots can then be utilized to optimally cover 
the target volume in minimum lime. The minimization of the 
duration of the treatment is an important factor for control- 
ling the cost of the procedure, especially when MR1 is used 
to guide and monitor the treatment. 16 If more control over 
the field pattern is required then an array wjih a larger num- 
ber of elements is required. Such an array could also be used 



Fig. 7. The input no^er requirement to reach the same temperature le*ei at 
the end of ultrasound pulsci for n anous cases. Case I ts lor the uniform 
excitation case, cast 2-5 arc corresponding to the amplitude and phase 
settings of cast t-iv in Table IV. 



to compensate for the effect of tissue nonuniforms es which 
may distort the focal pattern. 11 !: 
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1 . An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 

5 

A. a plurality of transducers, 

B. an excitation source for driving the plurality of transducers to generate and transmit 
ultrasound through the skull to induce cavitation at least at a selected region of the brain, 

10 

C. the excitation source driving at least selected transducers at differing phases with respect 
to one another. 

2. An. apparatus according to claim 1, wherein 

15 

A. the excitation source drives each of the selected transducers at a phase that compensates 
for a phase shift effected by the skull in the ultrasound generated and transmitted to the 
selected region by that transducer, 



20 B. so that the ultrasound generated by the selected transducers arrive substantially in phase 
with one another at the selected region. 



3 . An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
25 transducers arrive at phases within 90° of one another at the selected region. 
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4. An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
transducers arrive at phases within 45° of one another at the selected region. 

5 5. An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
transducers arrive at phases within 20° of one another at the selected region. 

6. An apparatus according to claim 1, wherein the excitation source drives the plurality of 
10 transducers to deliver ultrasound to the selected region at a frequency ranging from 0.01 

MHz to 10 MHz. 

7. An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0. 1 

15 MHz to 2 MHz. 

8. An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

20 

9. An apparatus according to claim 8, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

10. An apparatus according to claim 9, wherein the excitation source drives the plurality of 
25 transducers to deliver ultrasound to the selected region with burst mode operation, where 

burst mode repetition varies from 0.01 Hz to 1 MHz. 

11. An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 
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an ultrasound transducer, 

an excitation source for driving the transducer to generate and transmit ultrasound 
through the skull to induce cavitation at least at a selected region of the brain. 

An apparatus according to claim 11, wherein the excitation source drives the transducer 
to deliver ultrasound to the selected region at a frequency ranging from 20 kHz to 10 
MHz. 

An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.02 
MHz to 10 MHz. 

An apparatus according to claim 13, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

An apparatus according to claim 13, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

An apparatus according to claim 15, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising the steps of 



47 



WO 98/07373 



PCT/US97/14760 



placing a plurality of transducers in the vicinity of an exterior surface of the patient's 
skull, 

driving the plurality of transducers to generate and transmit ultrasound through the skull 
to induce cavitation at least at a selected region of the brain, 

the driving step including driving at least selected transducers at differing phases with 
respect to one another. 

A method according to claim 17, wherein step (C) includes driving each of the selected 
transducers at a phase that compensates for a phase shift effected by the skull in the 
ultrasound generated and transmitted to the selected region by that transducer, such that 
the ultrasound generated by the selected transducers arrive substantially in phase with 
one another at the selected region. 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 90° of one another at the selected region. 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 45° of one another at the selected region. 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at' a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 20° of one another at the selected region. 
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22. A method according to claim 17, wherein step (B) includes driving the plurality of 

transducers to deliver ultrasound to the selected region at a frequency ranging from 0.02 
MHz to 10 MHz. 

5 23 . A method according to claim 17, wherein step (B) includes driving the plurality of 

transducers to deliver ultrasound to the selected region at a frequency ranging from 0.1 
MHz to 2 MHz. 

24. A method according to claim 17, wherein step (B) includes driving the plurality of 

10 transducers to deliver ultrasound to the selected region at a sonication duration ranging 

from 100 nanoseconds to 30 minutes. 

25. A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

15 

26. A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

20 27. A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 

A. placing an ultrasound transducer in a vicinity of an exterior surface of the patient's skull, 

25 B. driving the transducer to generate and transmit ultrasound through the skull to induce 
cavitation at least at a selected region of the brain. 
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28. A method claim 27, wherein step (B) includes driving the transducer to deliver 
ultrasound to the selected region at a frequency ranging from 0.02 MHz to 10 MHz. 

29. A method according to claim 28, wherein step (B) includes driving the plurality of 

5 transducers to deliver ultrasound to the selected region at a frequency ranging from 0. 1 

MHz to 2 MHz. 

30. A method according to claim 29, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 

10 from 100 nanoseconds to 30 minutes. 

31. A method according to claim 30, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

15 32. A method according to claim 3 1, wherein step (B) includes driving the plurality of 

transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

33. An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
20 apparatus comprising 

A a plurality of transducers arranged in a two-dimensional array, and 

B. an excitation source for driving the plurality of transducers to generate and transmit 
25 ultrasound through the skull to induce cavitation at least at a selected region of the brain. 

34. An apparatus according to claim 33, wherein the plurality of transducers are arranged in 
an array having any of a substantially circular and a substantially annular cross-section. 
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35. An apparatus according to claim 34, wherein the plurality of transducers are arranged in 
an array having a substantially circular cross-section. 

5 36. An apparatus according to claim 34, wherein the plurality of transducers are at least one 
of (i) mounted in, and (ii) separated from one another by, a damping agent. 

37. An apparatus according to claim 36, wherein the damping agent is any of a natural or 
synthetic rubber. 

10 

38. An apparatus according to claim 34, wherein the excitation source drives at least selected 
transducers at differing phases with respect to one another. 

39. An apparatus according to claim 1, wherein excitation source drives each of the selected 
15 transducers at a phase that compensates for a phase shift effected by the skull in the 

ultrasound generated and transmitted to the selected region by that transducer, so that the 
ultrasound generated by the selected transducers arrive substantially in phase with one 
another at the selected region. 

20 40. A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising the steps of 

A. placing a plurality of transducers in the vicinity of an exterior surface of the patient's 
skull, 

25 

B. driving each of at least selected transducers generate and transmit ultrasound through the 
skull, and determining a phase shift effected by the skull in ultrasound generated and 
transmitted by each such transducer, 
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driving the plurality of transducers, together, to generate and transmit ultrasound through 
the skull to induce cavitation at least at a selected region of the brain, 

the driving step including driving at least selected the transducers at phases determined in 
accord with step (B) so that the ultrasound generated by at least the selected transducers 
arrive substantially in phase with one another at the selected region. 

A method according to claim 40, wherein the selected region ranges from 1 mm 3 - 1 cm 3 
in volume. 
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